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ABSTRACT 

Spitzer Space Telescope IRAC 3 — 8 /iin and AKARI IRC 2 — 4 /^m photometry are reported for ten white dwarfs with 
photospheric heavy elements; nine relatively cool stars with photospheric calcium, and one hotter star with a peculiar 
high carbon abundance. A substantial infrared excess is detected at HE 2221 — 1630, while modest excess emissions are 
identified at HE 0106— 3253 and HE 0307+0746, implying these latter two stars have relatively narrow (Ar < 0.1 i?0) 
rings of circumstellar dust. A likely 7.9 /im excess is found at PC 1225—079 and may represent, together with G166-58, 
a sub-class of dust ring with a large inner hole. The existence of attenuated disks at white dwarfs substantiates the 
connection between their photospheric heavy elements and the accretion of disrupted minor planets, indicating many 
polluted white dwarfs may harbor orbiting dust, even those lacking an obvious infrared excess. 

Subject headings: circumstellar matter — minor planets, asteroids — planetary systems - stars: 
abundances — white dwarfs 
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1. INTRODUCTION 

The large-scale composition of extrasolar minor planets 
can be indirectly m easured using white dwar fs as astro- 
physical detectors ( [Jura et al. 2009b; Zucker man et al.l 
[200 7). Rings of warm dust (and in some cases 
gaseous debris) revealed at over one dozen white dwarfs 
^arihi et al. 2009) are situated within the Roche limits 
([Gansicke et al.i i200Q) of their respective stellar hosts, 
and likely originated via th e tidal disr uption of minor 
(or possibly major) planets puralf200 3^. These closely 
orbiting disks rain metals onto the stellar photosphere, 
and contaminate an otherwise pristine hydrogen or he- 
lium atmospheric composition. For any particular cool 
white dwarf, the timescales for individual heavy elements 
to sink below the high gravity photosphere differ by less 
than a factor of three, yet are alwa ys orders of magni- 
tude shorter than the cooling age (Koester 2009a). With 
an appropriate treatment of the accretion history, photo- 
spheric heavy element abundances can be used as a mea- 
sure of the composition of the tidally destr oyed, pollu ting 
parent body or bodies ()Klein et al.ll2010l : lJurail2008D . 

The spectacularly debris-polluted white dwarf GD 362 
is an excellent example of this potential. Its circum- 
stellar disk re-emits at least 3% of the incident stellar 
radiation; one-third of thi s is carried by a s trong 10 ^m 
silicate emission feature (jJura eFaDIlQOTS). The star 
itself is polluted by at least 15 elements heavier than 
hydrogen and helium, in an arrangement that is rich in 
refractory elements and deficient in volatiles; a pattern 
that broadly mimics the inner Solar system, and compa- 
rable with the outer com position of the Earth and Moon 
(jZuckerman et aLll2007f l. A recent analysis of the un- 
usual atmospheric mix of hydrogen and helium, together 
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with X-ray constraints on the current hydrogen accretion 
rate at GD 362, suggests that if a single parent body gave 
rise to both the currently observed disk and the panoply 
of heavy elements in its atmosphere, then its mass would 
be larger than Ca llisto, and possibly larger than Mars 
pura et al.ll2009b[ ). 

Observations of a large number and variety of metal- 
polluted white dwarfs are necessary to better understand 
and constrain the connection between rocky extrasolar 
parent bodies, the circumstellar dust and gas, and the 
subsequ ent photospheric polluti ons they create in white 
dwarfs (jZuckerman et aLl I2003D . The lifetimes of the 
disks at white dwarfs are poorly constrained at present, 
though there is some indi cation their lifetimes do not 
signifi cantly exceed 10^ yr (jFarihi et al.ll2009l: IKilic et ahl 
120081 ). Another outstanding issue is the nature of the 
large metal abundances in stars lacking excess infrared 
emission; tenuous, gaseous, or previously-accreted disks 
are strong possibilities. While the evolution of optically 
thick dust at white dwarfs should be dominated by vis- 
cous dissipation-spreading, to date only G166-58 stands 
out as a possible indicator of how competing mechanisms, 
such as additional imp acts or gas drag at the inner disk 
edge, may play a role (|Jurall2C)08l: IFarihi et al.ll2008D . 

Spitzer IRAC and AKARI IRC photometric imaging 
observations are presented for nine metal-polluted white 
dwarfs chosen for their high calcium abundances and 
relatively short cooling ages. The observations are pre- 
sented in §2, the results in §3, and the implications for 
the minor planet accretion hypothesis and future dust 
searches are discussed in §4. 

2. OBSERVATIONS AND DATA 

The target sample of stars was chosen based on the 
prior result that nine of 14 white dwarfs with inferred 
metal accretion rates above 3 x lO^gs"^ and cooling 
ages less than 1.0 Gyr have infrared excess (jFarihi et al.l 
12001 . AU nine metal-contaminated white dwarfs in 
Table [T] satisfy this criteria; five DAZ stars with high 
rates of ongoing metal accretion, plus three DBZ stars 
with high calcium abundances. Table [2] lists the rel- 
evant parameters for the DAZ and DBZ stars. Also 
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Table 1 

Spitzer and AKARI White Dwarf Targets 



WD 


Name 


SpT 


V 

(mag) 


fCa/H(e)l 


Telescope 


Refs 


0047+190=^ 


HS 0047+1903 


DAZ 


16.1 


-6.1 


Spitzer 


1 


0106-328'' 


HE 0106-3253 


DAZ 


15.5 


-5.8 


Spitzer 


1 


0307+077=' 


HS 0307+0746 


DAZ 


16.4 


-7.1 


Spitzer 


1 


0842+231=' 


Ton 345 


DBZ 


15.9 


-6.9 


AKARI 


2 


1011+570 


GD 303 


DBZ 


14.6 


-7.8 


Spitzer /AKARI 


3 


1225-079 


PG 


DZAB 


14.8 


-7.9 


Spitzer 


3,4 


1542+182 


GD190 


DBQ 


14.7 




Spitzer 


5 


1709+230 


GD 205 


DBAZ 


14.9 


-8.0 


Spitzer / AKARI 


4,6 


2221-165" 


HE 2221-1630 


DAZ 


16.1 


-7.2 


Spitzer 


1 


2229+235" 


HS 2229+2335 


DAZ 


15.9 


-5.9 


Spitzer 


1 



Re fe rences. — 1 ) IKoester fc Wilken' '200E 2) 'Giinsicke et al. 200g 3) IWolff et al.l 
[200^ 4) IKoester et all l2005 5) Petitclerc et al. 2005 6) Voss et al. 2007; 
" The WD numbers for these stars are unofficial designations, but correctly reflect the 
conventional use of epoch B1950 coordinates. Abundances are expressed as [X/Y] = 
log [n(X)/n(Y)]. 



included as a target is the 22 000 K DBQ white dwarf 
GD 190 . with an anomalous atmospheric carbon abun- 
dance (|Petitclerc et al.l 120051 ). Ton 345 was included 
based on the presence of circumstellar, gaseous metal 
emission (jGansickc et al. 2008). 

2.1. Spitzer IRAC Observations 

Mid-infrared imaging observations of the white dwarf 
targets were obtaine d with the Spitzer Space Telescope 
(jWerner et al.ll20M ) during Cycle 5 us ing the Infrared 
Array Camera (IRAC: lFazio et al]l2004[ ) at all four wave- 
lengths, namely 3.6, 4.5, 5.7, and 7.9 ^m. The total in- 
tegration time in each channel was 600 s, where the ob- 
servations consisted of 20 frames taken in the cycling 
(medium) dither pattern with 30 s individual exposures. 
All images were analyzed, including photometry and up- 
per limits, as in Farihi et al. (2009) using both 1'.'2 and 
0'.'6pixel~^ mosaics created using MOPEX. 

In cases where the flux of a neighboring source was 
a potential contaminant of the white dwarf photometry, 
steps were taken to minimize or remove any such ex- 
ternal contributions, including small aperture radii and 
point-spread function (PSF) fitting with daophot. To this 
end, the newly available (beginning with pipeline ver- 
sion S18.5) 0'.'6pixel~^ mosaics were useful to better con- 
strain and eliminate sources of photometric confusion. 
Whil e the IRAC Data Handbook (Spitzer Sci ence Centeii 
|2006( ) provides information necessary to perform pho- 
tometry on point sources in the native l'.'2pixel~^ scale 
images, the provided aperture corrections do not take 
advantage of the available smaller pixel scale. A com- 
parison of identical observation sets, processed into both 
1'.'2 and 0'.'6pixel^^ mosaics, reveals a modest but clear 
increase (around 20% by full width at half maximum) 
in spatial resolution. This information gain at the 
smaller pixel scale is most evident at the two short- 
est wavelengths channels, where the diffraction lim- 
ited PSF of the tele scope is significantly under-sampled 
(jSpitzer Science C"en tcr 2006). 

2.2. AKARI IRC Observations 

Observations with AKARI ijMurakami et all I2007D 
were executed u sing the Infrared Camera (IRC; 
lOnaka et al.) [Ml in the N2, N3 and N4 fihers (2.34, 



Table 2 

Spitzer IRAC Fluxes and Upper Limits 



WD 














(MJy) 


(MJy) 




0047+190 


73 ±4 


47 + 3 


28 + 4 


16 ± 5 


0106-328 


151 + 8 


106 + 6 


74 + 5 


55 ± 5 


0307+077 


103 ± 5 


72 + 4 


54 + 5 


45 ±6 


1011+570 


220 ± 11 


137 + 7 


86 + 6 


55 ±5 


1225-079 


328 ± 16 


224 ± 11 


139 + 8 


106 + 8 


1542+182 


194 ± 10 


121 + 6 


75 + 6 


38 ±5 


1709+230 


180 + 9 


107 + 6 


63 + 5 


32 ±5 


2221-165 


204 ± 10 


172 + 9 


157 + 9 


119 + 8 


2229+235 


61 + 6 


37 + 4 


20 + 5 


18" 


" 3(T upper 


limit. 









Table 3 

AKARI IRC Fluxes 



WD 




Fz.2iJ.ui 


F4.z1.1m 




(MJy) 


(MJy) 


(MJy) 


0842+231 


293 ± 22 


271 ± 23 


218 ± 18 


1011+570 


586 ± 46 


289 ± 22 


146 ± 12 


1709+230 


418 + 32 


237 ± 18 


111 ± 11 



3.19, and 4.33 /im) for three white dwarfs listed in Ta- 
ble [TJ These stars were imaged in 2008 November and 
2009 March- April using the Astronomical Observational 
Template IRCZ3, which was designed for imaging with 
the N2, N3, and N4 filters in a pointed observation. The 
basic data reduction, including dark subtraction, linear- 
ity correction, distortion correction, flat fielding, image 
stacking, and absolute position determination were per- 
formed by the standard IRC Imaging Data Reduction 
Pipeline for Phase 3 data (version 080924P3). The IRAF 
package daophot was used for photometry, with an aper- 
ture radius of 5 pixels (7'.'3) and corrected to the standard 
10 pixel radius using aperture corrections derived from 
bright sources within individual images. Flux calibra- 
tion was achieved using conversion factors provided by 
the IRC team^ 

^ http:/ /www. ir. isas.jaxa.jp/ASTRO- 

F/Observation/DataReduction/IRC/ConversionFactor__090824.html 
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Figure 1. SEDs of HS 0047+1903, HE 0106-3253,HS 0307+0746, and Ton 345 fitted with DA and DB white dwarf models (iKoesteil2009y). 
Filled circles represent the photometry (Melis et al. 2010; Brinkworth et al. 2010; Farihi 2009; Abazaiian et al. 2009; Skrutskie et al. 200l; 
[Uopenhagen University Obser vatory., 2006 : .DENIS Consortium. 2005, ; Martin et al.,,,2005.) while open circles are fluxes based on spectroscop- 
ically derived colors llK'TTic'er^l.ll2006i) , and hence not independent measurements. 



3. ANALYSIS AND RESULTS 

Tables list the flux determinations and established 
upper limits for the science targets. Figures [T]-[6] plot 
the measured mid-infrared fluxes toge ther with avail- 
able near-infrared and op tical da ta (Abazaiian ct al.! 
20pl Farihi 2009; Skrutskie et all '2006; Martin ct al. 
2005; McCook & Sion 2008; Landoh & Uomoto 2007; 
Kilic et al.i i2006i: Copenhagen University Observatoryi 



2006t ^ I Space Te lescope Scie nce Institute! 



Zacharias et al 



2005i: iDENIS Consortium! 



2006 



2005 



Monet ct all 2003). For a few of the metal-rich 
target stars, !Kilic et all (j2006[ ) has published spectro- 
scopically derived J — H and H — K colors, such that the 
H- and fiT-band fluxes of these stars are anchored to the 
2MASS J-band flux, and hinge upon its accuracy. IRTF 
JHK photometry of HE 0106-3252 and PC 1225-079 
are taken from iFarihil ([2009|). Where availa ble, Galaxy 
Evolution Explorer {GALEX; !Martin et all 120051 far- 
and near-ultraviolet photometric fluxes are plotted to 
assist in establishing the appropriate brightness level 
of the white dwarf. However, because the interstellar 
extinction to these stars is not known, the GALEX 
fluxes were used with some caution. 

Five of the targ et white dwarfs co me from the Ham- 
burg Schmidt (HS; Hagen et al.!!l995D a nd the Hamburg 
ESO Schmidt (HE; 



Wisotzki et al.!!l996l ) quasar surveys. 



The coordinates in SIMBAD for these five white dwarfs 
have poor precision , and are o ften missing from the white 
dwarf catalog McCook fc Sion (1999, 2008). The correct 
positions for these white dwarfs were provided by the 
SPY team (R. Napiwotzki 2007, private communication), 
and their coordinates in the IRAC images are given in 
Tabled) The HE and HS white dwarfs tend to suffer from 
a lack of accurate ground-based photometry in the liter- 
ature. A typical star has entries in a few photographic 
catalogs (30% errors), variable quality r'-band CCD pho- 
tometry from the Carlsberg Meridian Ca talog (CMC; 
!Copenhagen University Observatory! !2006! ). and a DE- 
NIS /-band measurement if in the southern hemisphere. 
For HS 0047+0093, HS 0307-^0746, HE 2221-1630, and 
HS 2229+2335, the only near-infrared photometry is the 
2MASS J-band value, as the H- and i^^-band data are 
subst antially less reliable ow ing to the faintness of these 
stars (jSkrutskie et al.ll2006D . 

All available photometric data for each white dwarf 
were fitted by a pure hydrogen or pure helium white 
dwarf model with log [g (cm s~^ )] = 8.0 and effective tem- 
perature to the nearest 500 K (|KoesteHl2009b) . Because 
the surface gravity and the atmospheric composition are 
fixed, and because the effective temperature resolution 
is limited, the models shown in the figures do not neces- 
sarily describe accurately or strongly constrain the effec- 
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Figure 2. SEDs of GD 303, PG 1225-079, GD 190, and GD 205 fitted DB white dwarf mo dels l|Koestegi2009bt). Filled circles represent 
tlie photometry {F arihi 2009; Abazaiian et al. 2009; Landolt k. Uomoto 2007; Skrutskic ct 'aIll200g : IDENIS Consortiumll2005l : IMartin et al.l 
[20051; [McCook fc Sioail999l) . The apparent 2MASS /Cs-band excess at GD 190 is spurious. 



Table 4 

IRAC Coordinates for HE and HS White 
Dwarfs 



Star 


a 


<5 




HS 0047+1903 


00''50™ 


12.47" 


+19°19 


48" 2 


HE 0106-3253 


0l''08'" 


36.04= 


-32°37 


43" 5 


HS 0307+0746 


03'' 10™ 


09.12= 


+07°57 


31'.' 5 


HE 2221-1630 


22h24m 


17.44= 


-16°15 


47" 7 


HS 2229+2335 


22''31™ 


45.49= 


+23°51 


24" 1 


Note. — As 


measured on the 


IRAC 


array 



based on image header astrometry; the epochs 
all correspond to Cycle 5 dates between 2008.9 
and 2009.3 (equinox 2000). 

tive temperature of the stars. Rather, the atmospheric 
model fits predict the expected level of mid-infrared pho- 
tospheric emission, and there by constr ain the presence or 
absence of infrared excess (Farihi 2009). Relative to the 
Spitzer data, the AKARI images were found to be shal- 
low, with patchy sky background, and highly asymmetric 
stellar images. For this reason, it is likely the errors used 
here are underestimated, and these data were given little 
weight in the final analysis. 

3.1. Stars With Infrared Excess 

HE 0106-3253. The IRAC images of this star reveal 
another source separated by 4'.'4 at position angle 299°, 



Table 5 

Flux Determinations! ^j. HE 0106-3253 



Method^ 




.F4.5nm 












(MJy) 




a 


154.4+ 10.9 


106.6 + 8.3 


75.3 + 7.6 


56.1 + 7.7 


b 


158.0+ 11.2 


108.0 + 8.4 


77.4 + 7.8 


56.1 + 7.7 


c 


147.1 ± 10.4 


106.4 + 7.6 


72.7 + 5.9 


54.0 + 5.5 



! Adopted fluxes are the weighted average of methods a and c. 
^ Methods: a) PSF fitting of both white dwarf and neighboring 
source; b) PSF fitting of white dwarf only; c) Relative aperture 
photometry at r = 1'.'2 using nearby star (§3.1). 



as measured by daophot using the 0'.'6 pixel""'^ images 
(Figure |3]). The neighboring source appears to be a back- 
ground galaxy based on its half maximum intensity at 3.6 
and 4.5 /xm, measured in a direction orthogonal to its po- 
sition angle relative to the white dwarf. Two independent 
methods were used to measure the flux of the white dwarf 
without contamination from the nearby galaxy. First, 
daophot was employed to fit both the white dwarf and 
the nearby galaxy (thus treating it as a point source) si- 
multaneously. This was done with two independent sets 
of PSF stars, one derived from archival IRAC images 
taken for the Galactic component of the First Look Sur- 
vey, and the other created using a few to several of bright 
point sources within the science target images. Addition- 
ally, daophot was run to fit the white dwarf only, with 
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a) 3.6 /j.m 
WD 

m ' 


b) 

• 


c) 

* 

V 


d) 

WD 1^ 

* 



Figure 3. IRAC 3.6 mosaics of HE 0106-3253: a) white dwarf 
and nearby source; b) both sources fitted and subtracted with 
daophot; c) white dwarf removed by daophot; d) expected white 
dwarf photospheric flux removed by daophot, revealing a point-hke 
excess. Each panel is 1' on a side where upwards corresponds to 
position angle 36.0°. 



a) 5,7 fj.m 


b) 




WD 






# 




• 


c) 


d) 




* 




WD ' ^ 


♦ 







Figure 4. Same as Figure 5, but for 5.7 fim. Again, excess flux is 
revealed after removal of the expected white dwarf photosphere. 

no information provided on the neighboring source. It 
is noteworthy that the task returned similar resuhs for 
these two approaches, within 3% at aU four IRAC wave- 
lengths. 

Second, relative aperture photometry was performed 
on HE 0106— 3253 and a moderately bright, nearby star 
2MASS J01083497-3237486; this star can be seen in Fig- 
ures [3] and |31 33'.'6 distant from the white dwarf. Robust 
centroids were determined for each star using both gaus- 
sian profile analysis and daophot PSF fitting; the cen- 
troids in all four channels were averaged for each star 
and found to agree to within 1 standard deviation of 
0.10 pixels. The centroid for the comparison star was 
taken to be fixed while photometry for the white dwarf 
was performed at a total of nine positions; at the nominal 
centroid plus at eight positions along a circle of radius 
0.10 pixels about this center, spaced evenly at angles 45 
to 360° in the xy image plane. Photometry was executed 
in half-pixel steps from 1 to 4 pixels. The same relative 
photometric analysis was also performed for two ran- 



domly selected point sources (called 'starl' and 'star2') 
within the HE 0106-3253 IRAC field of view. Lastly, the 
2MASS comparison star was input to daophot as a refer- 
ence PSF, and the white dwarf flux was re-measured as 
above. 

Figure [S] plots the results of the relative flux measure- 
ments between the white dwarf and 2MASS comparison 
star. Each datum is the average and 3(7 error among the 
nine positions where the white dwarf flux was measured. 
Also shown in the plot are the fluxes as determined by 
daophot, averaged over the 3 sets of PSF stars, and the 
expected white dwarf photospheric fluxes as determined 
by the model shown in Figure [TJ Two things are ap- 
parent from the figure and analysis: 1) HE 0106— 3253 is 
photometrically extended beyond r = 2 pixels, and possi- 
bly beyond 1.5 pixels; 2) the daophot flux measurements 
are in good agreement with the fluxes determined by rel- 
ative photometry with the 2MASS comparison star at 
r — 2 pixels. 

The fluxes determined by these methods are listed in 
Table [5l and a series of representative images are shown 
in Figures [3] and |31 Both methods yield IRAC fluxes 
in excess of the predicted photospheric emission of the 
white dwarf based on the model plotted in Figure [TJ 
the weighted average of these is plotted in the Figure 
and listed in Table [21 The excess at each IRAC wave- 
length lies in the range 3 — 5cr using the errors derived 
via the individual measurements, or within 5 — 6(T using 
the weighted average errors. The nearby galaxy can- 
not account for the results of the PSF fitting photome- 
try where extracted sources must conform to a relatively 
static, model point source. The last panel of Figures [3] 
and [4] displays the PSF fit and removal of the white dwarf 
at the expected photospheric level; a point-like residual 
(excess) is seen at the location of the white dwarf. 

To assess the potential flux contamination by the 
galaxy, a radial profile and vector analysis of the white 
dwarf fltix was performed. Both ID gaussian profile fit- 
ting and a raw vector cut (each of 3 pixel width) were 
executed across the white dwarf in the direction towards 
and away from the galaxy, and the fluxes on either side 
were compared. From these analyses, it is found that 
at r = 2 pixels from the white dwarf centroid, the side 
nearest the galaxy contains at most 15% more flux than 
the side opposite at 3.6 /im, and no more than 10% at 
the three longer wavelengths. Assuming a worst case 
scenario where one-half of the photometric aperture is 
biased in this manner - certain not to be the case - then 
the galaxy contributes a maximum of 5% of the measured 
flux at all wavelengths except 3.6 fiui, where it is limited 
to 7.5%. The excess flux (measured— photosphere) rela- 
tive to the measured flux is 27, 32, 39, and 55% at 3.6, 
4.5, 5.7, and 7.9 /im, and hence the potential influence of 
the galaxy cannot account this excess. From this analy- 
sis, it appears that HE 0106—3253 has an infrared excess. 

HS 0307+0746. Figure [1] reveals an infrared excess 
around this star based on the plotted 10 500 K DA model 
and the short wavelength photometry. The HK data 
are anchored to the 2MASS J-band measurement of 
16.39±0.14mag, where the J — H and H — j'C colors were 
derived spectroscopically bv lKilic et al.l (|2006[ ). The only 
optical, photoelectric photometry comes from the Carls- 
berg Meridian Catalog (CMC) where two measurements 
are given as r' = 16.48 ± 0.08 mag. .Koester fc WilkenI 
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Figure 5. IRAC fluxes for HE 0106-3253 derived by photometry relative to the nearby comparison star 2MASS J01083497-3237486 as 
a function of aperture radius. The data points and error bars represent the average and 3cr deviations for an ensemble of nine white dwarf 
centroids. The dotted line marks the white dwarf flux as determined by daophot, while the dashed line shows the expected photospheric 
flux. 



()2006D give T^e = 10 200 K and log g = 7.9 for the white 
dwarf, similar to the pfotted model that agrees well with 
the GALEX iai- and near-ultraviolet fluxes. If the model 
used here is sufficiently accurate, the putative excess at 
HS 0307+0746 appears at the 4 - 6ct level at aU IRAC 
wavelengths. 

Ton 345. The Spitzer data for this star is analyzed 
elsewhere, but the AKARI IRC fluxes reported here 
are original and appear broadly co nsistent with the 
near-infrared ami Spitz er IRAC fluxes (jMelis et al.ll201Cil : 
iBrinkworth et al.ll2010D . 

PG 1225-079 and GD 303. The 7.9 ^fm IRAC flux at 
PG 1225-079 indicates a Aa excess (Figure [2|) while ah 
other photometric data are consistent with photospheric 
emission. If real, this star joins G166-58 as a metal-rich 
white dwarf with a long wavelength IRAC excess but 
no short wavelength excess; these two objects may form 
a class of object whose orbiting disks have large inner 
holes. GD 303 may have a marginal photometric excess 
at the 2cr level (Figure [5]). 

HE2221-1630. This metal-polluted white dwarf 
shows the clearest mid-infrared excess (Figu re of all 
the science targets. iKoester fc WilkenI ()2006[ ) give Teg — 
10 100 K for this star. 

3.2. Notes on Stars Without Infrared Excess 



GD 190. The 2MASS A% = 14.96 ± 0.13 mag measure- 
ment (Figure [5]) for this helium- and carbon-rich white 
dwarf appears erroneous at the 2.5cr level. 

GD 205. The fit to the spectral energy distribution 
(SED) shown in Figure [5] employs a hi gher T^s tha n that 
found by previous stud i es; 1 8 500 K (jVoss et al.i r2007'). 
19 700 K C Koester et al|[2005 :). However, the 20 500 K 
model better matches all the photometric data, including 
the GALEX Huyies. 

2229+235. While the native pixel scale IRAC images 
of this star reveal a single source (single Airy disk) with 
an elongated core, the 0'.'6 pixel"'^ mosaics spatially re- 
solve a neighboring source at 2'.'0 and position angle 336°, 
as measured by daophot. The 2MASS J- and H-hand im- 
ages of the white dwarf may show a slight elongation in 
the expected direction, but the S/N in these images is 
too low to be certain (especially at _ff-band) . Photomet- 
ric deconvolution was performed with daophot, yielding 
good results at 3.6 and 4.5 fim; errors of 6% and 8% 
respectively for the white dwarf. While the neighboring 
source is likely extragalactic in nature, there is no evi- 
dence of elongation in the IRAC images, and there are 
no significant residuals in the PSF subtracted images of 
both sources at either of the shorter wavelengths. 

Unfortunately, reliable optical or near- 
infrared photometry of HS 2229+2335 is lack- 
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; HE 2221-1630 



T = 10000 K 



0.01 I 

0.1 



1.0 10.0 



HE 0106-3253 



HS 2229+3335 



T = 18500 K 



0.1 1.0 10.0 100.0 

A (^m) 

Figure 6. SEDs of HE 2221-1630 and HS 2229+2335 fitted with 
DA wiiite dwarf models (Koester 2009b). Filled circles repre- 
sent the photometry ( Copenhagen University Observatory '200^; 
Skrutskic et al. 2006; Space Telescope Science Institute 200B; 
DENIS Consortium 2005; Zacharias et al. 2005; Martin et al. 2001; 
Wonet et al. 2003), while open circles arc fluxes based on spectro- 
scopically derived colors iKilic et al.;,.2006 ) , and hence not indepen- 
dent measurements. The downward arrow is a 3a upper limit. 

ing; the best data are r' = 16. 18 ± 0.25 mag 
(jCopenhagen University Obs ervatorV '2006') , and 
J = 16.16 ±0.09 mag (jSkrutskie et al. 2006). The 
flux plotted at F-band is the average and stan- 
dard deviation of five photographic plate magni- 
tudes from NOMAD, GSC2.3 and USNO-Bl.O 
(I Space Telescope Science Institutell2006l : [Zacharias et al.l 
120051 iMonet et al.l \2m^ . iFinterstellar extinction 
towards the white dwarf is negligible, the GALEX and 
IRAC data agree w i th an 18 500 K, log g = 8.0 DA model 
iKoester fc Wilke"nl ()2006[ ). Near- infrared photometry 
would better constrain the SED of this white dwarf, but 
it does not appear to have a mid-infrared excess. 

3.3. Disk Models 

The infrared excesses at HE 0106— 3253, 
HS 0307-^0746, and HE 2221-16 30 w e re fit ted with 
the circumstellar disk models of lJural ()2003[ ). These 
disks are geometrically thin, flat rings of dust vertically 
optically thick at wavelengths up to 20 /im. Figure [7] 
shows the model fits, while Table |6] gives the model 
parameters. The models do a reasonable job of re- 
producing the disk fluxes. Both HE 0106— 3253 and 
HS 0307-1-0746 have a relatively high 3.6 /im flux relative 
to the model; this behavior is also seen for the inner 



HS 0307+0746 



4 6 
A {jim) 



HE 2221-1630 



A (^m) 

Figure 7. The simple flat ring model applied to the excess fluxes 
at HE 0106-3253, HE 0307+0746, and HE 2221-1630. 



disk region at PG 1457-086 (|Fari^etaL| [200|) . It is 
conceivable that some disk warping in the inner regions 
is responsible for the relatively high observed flux at 
short infrared wav elengths, similar to tha t seen and 
modeled for CD 56 iTura et al.ll2m9aL l2007al) . 



4. DISCUSSION 



4.1. 



The Existence of Narrow Dust Rings 

Two white dwarfs in this study - HE 0106— 3253 and 
HS 0307+0746 - appear to have narrow (Ar < 0.1 i?©) 
rings of dust similar to PG 1457—086, when compared 
to the en semble of white dw arfs with warm circumstel- 
lar disks (jFarihi et al.l 120091 ). In fact, their disk widths 
could be even smaller as the models presented here are 
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Table 6 

White Dwarf Flat Ring Model Parameters 



WD 




R/d 


T- 

^ inner 


Toutcr 


"^inncr 


'^outcr 


cos i 




(K) 


(10-12) 


(K) 


(K) 


(fi©) 






0106-328 


15 700 


3.4 


1400 


1100 


0.19 


0.27 


0.15 


0307+077 


10 200 


3.5 


1050 


850 


0.17 


0.22 


0.40 


2221-165 


10100 


4.3 


1200 


750 


0.14 


0.27 


0.50 



Note. — The first column is the input temperature of the central 
star, followed by the ratio of the stellar radius to its photometric 
distance. 



degenerate between disk width and inclination. For ex- 
ample, if the inclination of the disks at HE 0106— 3253 or 
PG 1457— 086 are actually close to face-on (z = 0°), the 
radial extent of their dust is only 0.01 Rq, roughly 1 i?® 
and 10 times smaller than the rings of Saturn, yet still 
over 10 times the diameter of all Solar System asteroids 
save Ceres. This finding has an important implication. 

The detection of narrow dust rings suggests the as- 
teroid accretion model may apply to additional, and 
potentially many, stars without a clear infrared ex- 
cess. Take the most modest excess yet detected, that 
of PG 1457-086: a face-on disk of width 0.01 i?© and 
producing a 15% (3ct) excess would become difficult to 
confirm above i — 50°, corresponding to an excess of 
only 10%. While such disks may appear slight when 
compared to the majority of infrared excesses discovered 
so far, they could easily contain sufficient mass to fuel 
the observed photospheric pollutions in DAZ and DBZ 
stars. At the high densities necessa ry to prevent coll i- 
sional dust erosion within the disk (jFarihi et al.ll2008D . 
the thin ring above could harbor over 10^^ g of matter 
within a disk height of 10 m, and readily supply heavy 
elements at lO^gs"^ for nearly 10^ yr. The discovery of 
three relatively narrow rings indicates there may be more 
subtle excesses awaiting detection. 

Table [7] lists a comparison of the infrared excesses at 
all known white dwarfs with dust. To calculate the frac- 
tional dust continuum luminosity, r = Lir/Lwd, the 
excess flux in the 2 — 6 /xm region was fltted by a single 
temperature blackbody, and its bolometric flux divided 
by the stellar flux (similar to those plotted in the SED 
figures). There are currently three white dwarfs with T ~ 
1000 K dust whose infrared excesses are less than 10^^ 
of the incident starlight: HE 0106-3253, WD 1041-^091 
(SDSS J104341.53+085558.2), and PG 1457-086. The 
stars G166-58 and PG 1225—079 owe their modest ex- 
cesses to a lack of T > 500 K dust, hence their low r 
are qualitatively different. It is noteworthy that all stars 
with T > 0.02 had their disks discovered via ground- 
based observations. 

A disk of radial width 0.01 Rq may result from the tidal 
disruption of an asteroid perturbed into an orbit with 
periastron near 0.2 i?Q, and dust spreading timescales 
suggest the material need not expand to a much greater 
radius. In the context of Spitzer searches for dust at 
polluted white dwarfs, highly attenuated rings of dust 
may persist at some white dwarfs that lack apparent 
infrared excess. Without supplementary ground-based, 
near-infrared photometry, the mild infrared excess at 
PG1457-08 6 would not be obvious (jFarihi et al.l[200l 
iFarihil 120091 ) . A combination of accurate optical and 



near-infrared photometry with appropriate atmospheric 
models may reveal excesses previously missed. Ideally, 
high signal-to-noise mid-infrared spectroscopy spanning 
the 10 /im silicate feature should be capable of detecting 
very modest dust emission, such as that seen at HD 69830 
(jLisse et al.l 120071) . Spectroscopy with the Mid- Infrared 
Instrument for JWST would be an ideal choice to reveal 
the most subtle infrared excesses. 

4.2. A Class of Disks with Enlarged, Dust-Free Inner 
Regions ? 

Among white dwarfs with circumstellar dust, G166-58 
was the first star observed to h ave an beginning long 
ward of 4 /im (iFarihi et al.l 120081 ) . Figure |8] revisits the 
available photometry for this DAZ star and plots the 
IRAC data, together with IRTF JHK, SDSS ugriz, and 
GALEX near-ultraviolet photometry. These data are fit- 
ted with a 7500 K, log g = 8.0 DA model that reproduces 
the photometry well to 4.5 /im. The IRAC fluxes have 
changed marginally following analysis of the newly avail- 
able 0'.'6pixel~^ mosaics for this star. The excess can 
be reproduced by a 500 K blackbody, implying the inner 
dust disk edge cannot be significantly warmer. The flat 
ring model then predicts the innermost dust is located at 
0.28 i?0, and far from where silicate grains should rapidly 
sublimate. Therefore, the dust disk appears to have an 
enlarged, dust-free inner region. 

Interestingly, PG 1225—079 appears to have an excess 
only at 7.9 /im and may also represent a dust disk with 
a large inner 'hole'. An excess flux is difficult to con- 
flrm with the current data alone, yet the IRAC images 
are of good quality and free of neighboring sources. If 
the excess is confirmed, the dust-free, inner region at 
PG 1225—079 should be around three times larger than 
that inferred for G166-58, but still located well within 
the Roche limi t for kilometer s ize or larger solid (rocky 
or icy) bodies (|Davidssonlll999D . There are two distinct 
possibilities for such large dust-free zones: 1) the region 
is dominated by gaseous matter emitting inefficiently in 
the infrared or, 2) the region is depleted relative to the 
dust inferred from the observed infrared excess. 

The first possibility can arise if the inner dust disk 
region is impacted by a moderate size asteroid whose 
solid mass and angular momentum is dis sipated via col- 
lisions prior to vaporizing the entire disk p ural 120081 ) . In 
Saturn's rings, the observed radial dark lanes or spokes 
are the result of meteorite impacts, whose kinetic ener- 
gies are sufficient to produce fine dust, gas, and (impor- 
tantly) plasma, which interact d ynamically within this 
environment to obscure the rings (jCoertz fc Morfilll[T983l: 
[MorflU fc Goertz..l983i) . Similar impacts upon dust disks 
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Table 7 

Thermal Continuum Excess at White Dwarfs 



nth Dust 



WD 


Name 


Stellar Model 2 - 
TwD (K) 


- 6 /im Blackbody" 
TiR (K) 




0106-328 


HE 0106-3253 


16 000 


1400 


0.0008 


0146+187 


GD16 


11500 


1000 


0.0141 


0300-013 


GD 40 


15 000 


1200 


0.0033 


0408-041 


GD 56 


14 500 


1000 


0.0257 


0307+077 


HS 0307+0746 


10 500 


1200 


0.0028 


0842+231 


Ton 345 


18 500 


1300 


0.0048 




PG 


19 500 


1200 


0.0017 


1041+091 


SDSS 104341.53+085558.2 


18 500 


1500 


0.0008 


1116+026 


GD133 


12 000 


1000 


0.0047 


1150-153 


EC 11507-1519 


12 500 


900 


0.0216 


1226+110 


SDSS 122859.93+104032.9 


22 000 


1000 


0.0042 


1225-079: 


PG 


10 500 


300: 


0.0005: 


1455+298 


G166-58 


7500 


500 


0.0015 


1457-086 


PG 


20 000 


1800: 


0.0007 


1729+371 


GD 362 


10 500 


900 


0.0235 


2115-560 


LTT 8452 


9500 


900 


0.0092 


2221-165 


HE 2221-1630 


10100 


1000 


0.0076 


2326+049 


G29-38 


11500 


1000 


0.0297 


Note. — 


Measured infrared excess from thermal continuum emission between 2 and 6 ^m, as 



most stars lack spectroscopic data on any potential silicate emission in the 8 - 
^ This single temperature is a zeroth order approximation of the disk SED. 



12 fim region. 



at white dwarfs should release kinetic energies 1000 times 
higher due to the difference in gravitational potential; the 
model radii and sizes of dust rings at white dwarfs are 
very similar to those of Saturn. 

This model is attractive because it can account for 
disks with enlarged, dust-free inner regions or small outer 
radii (including narrow rings) as well as high accretion 
rate DAZ white dwarfs that lack infrared excess out to 
7.9 fim yet require a nearby, ongoing source of photo- 
spheric heavy element pollution. G166-58 has a relatively 
long diffusion timescale fo r a DAZ star, a round 2000 yr 
for the observed calcium ()Koesteiil2009al) . This calcu- 
lated lifetime does not guarantee accretion is ongoing, 
but the presence of orbiting dust makes it likely. 

Perhaps more important, the absence of detectable 
dust need not imply an absence of accretion and pho- 
tospheric pollution. While the disks at G166-58 and 
PG 1225—079 may be partly composed of dust par- 
ticles which are detected in the infrared, and partly 
of gaseous heavy elements which remain unseen, most 
metal-contaminate d white dwarfs lac k infrared dust sig- 
natures altogether (jFarihi et al.l l2009). In many of these 
cases accretion must be ongoing and can be accounted for 
by the presenc e of either t enuous dust rings or completely 
gaseous disks (|Jurall2008D : unseen bro wn dwarf compan- 
ions have been ruled out (jFarihi et al . 2009, 2008). 

The second possibility for dust disks with large inner 
holes is a relative dearth of dust (and gas) due to prior ac- 
cretion; i.e. a partly consumed disk. PG 1225—079 has a 
heavy element diffusion timesc ale typical of co oler DBZ 
stars, over 1 Myr for calcium (jKoesteii l2009al) . In this 
case the lifetime of metals implies that accretion need 
not currently proceed at a high rate, and the observed 
metals may be residuals from an older event. Although 
the presence of an infrared excess, and hence the possi- 
bility of orbiting dust, make an ongoing accretion sce- 
nario likely, it raises an important question about the 
disk depletion timescale at white dwarfs. A prediction 
of ceased accretion would be higher relative atmospheric 



abundances of lighter elements (jKoesteii l2009aj) . and is 
testable with good abundance data for PG 1225—079. 

In the appendix, a simple calculation demonstrates 
that radiation drag (the Poynting-Robertson effect) on 
particles is insufficient to drive the infall of disk material 
and account for the distribution of inferred h eavy element 
accretion rates at disk- polluted DAZ stars (iFarihi et alj 
[200l[j^IraeFaI1l2007bl:IKoester fc Wilkenll2006D . At the 
inner disk edge, viscous spreading due to collisions be- 
tween differentially rotating r ing an nuli and random mo- 
tions among solids (Espositolll993| ) must be boosted by 
the presence of gas (via sublimated or coUisionally va- 
porized dus t ); suc h enhancement must be significant. 
iFarihi et al.l ( 2008D have shown that grains within white 
dwarf disks will collide on a timescale at least an order 
of magnitude shorter than their Poynting-Robertson life- 
time. Therefore, high viscosity resulting from dust-gas 
interactions must be ultimately responsible for the im- 
plied accretion rates from observed disks. For a disk that 
is otherwise optically thick at all relevant wavelengths, 
only a tiny fraction of its total volume - only the inner 
edge ~ will be illuminated by the central star. There 
is no compelling reason why the remainder of the disk 
should not evo lve as a particu late ring, having a very low 
dust viscosity (jEspositolll993D , and be gradually whittled 
away via gas drag at the inner edge. 

4.3. The Emerging Picture 

Figures [9] and [10] plot the metal abundances versus ef- 
fective temperature and time-averaged metal accretion 
rate versus cooling age for 60 Spitzer observed metal- 
rich white dw arfs. All accr etion rates were re-calculated 
following Fari hi et al.ll2009l using recent results on metal 
diffusion timescales and convective envelope masses for 
both DA Z and DBZ white dwarfs, including mixed cases 
(jKoested r2009a: D. Koester 2009, private communica- 
tion); parameters for the newly observed Cycle 5 stars 
are listed in Table [S] Plotted as newly discovered disks 
are HE 0106-3253, HE 0307-^0746, HE 2221-1630, and 
tentatively PG 1225-079. To these are added three 
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0.01 I 
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Figure 8. SED of G166-58 fitted with a DA wliite dwarf 
model (Kocstcr 2009b;^ . Filled circles represent tlie p liotometry 
Pbazaiian et al..,2009l : IFarihi|[2009l: IMartiii et al|[2005f) . 



DAZ 
DBZ 

DAZ, disk 
DDZ, disk 



15 

T (10'' K) 



Figure 9. Measured calcium abundance versus effective tempera- 
ture for all 61 metal-rich white dwarfs observed with Spitzer IRAC. 
Abundances are expressed as [X/Y] = log [n(X)/ra(Y)]. 

whit e dwarfs orbited by bo th dust and gaseous met- 
als (jBrinkworth et aLI [2009t C. Brinkworth 2008, pri- 
vate communication), while PG 1632 +177 (observed by 
iFarihi et alll2008l and listed as DAZ in lFarihi et "alll2009l) 
has been re moved from this growin g sample, as it is not 
metal- lined (jZuckerman et aLll2003[ ). 

The new additions to the plots continue the trend 
that metal accretion rates dM/dt ^ 3 x lO^gs"^ are 
those that are rn ost likely to exhibit infrared excess 
(|Farihi et al.ll2009D : the sole exception is G166-58. Two 
of the disks discovered in this work orbit relatively cool 
white dwarfs with cooling ages of 0.73 and 0.94 Gyr. 
These are somewhat older than the previous average 
cooling age of w hite dwarfs with disks as found in 
IFarihi et al.l (pOOl ). and may indicate relative longevity 
for orbiting dust. 

Viscous spreading among solids may allow optically 
thick dust rings at white dwarfs to persist for at least 
lOMyr timescales, possibly orders of magnitude longer, 
analogous to planetary rings. The inner radius of Sat- 
urn's rings represents the furthest inne r extent of its 
sprea ding over a period up to a few Gyr ()Esposito et al.l 
I2008D . For dust rings at white dwarfs, gas drag should de- 
crease the viscous timescale by orders of magnitude, but 
the gas content in most white dwarf disks is presently 



9.0 
log [t„ 



(yr)] 



Figure 10. Time-averaged dust accretion rates versus cooling age 
for all 61 metal-rich white dwarfs observed with Spitzer IRAC. The 
dotted line represents an accretion rate of 3 X 10** gs^^. Table [8l 
lists the {x, y) coordinates for the new dust disk identifications. 

unconstrained. At the inner edge of a typical dust disk, 
rapid sublimation of dust will produce significant g as, as 
may grain-grain collisions. Such viscosity enhancements 
are necessary to reproduce the dust accretion rates in- 
ferred for stars with short photospheric metal lifetimes. 

Dust rings at white dwarfs are unlikely to be as simple 
as modeled here; calcium triplet emission lines at three 
metal-polluted stars can be approximately reproduced by 
orbiting disks, but there are asymmetries which indicate 
deviations from the model (Gansicke et al. 2008). Yet 
it is noteworthy that the simple model works rather well 
given only three free parameters, two of which - the inner 
and outer radius - are physically well-motivated. 

Significant orbital energy must be shed to transform an 
asteroid on a perturbed, eccentric orbit into a disk con- 
tained entirely within the Roche limit of a white dwarf. 
The primordial energy may dissipate partially via colli- 
sions which initially widen the disk (or perhaps not in 
the case of narrow rings), yet an efficient way to discard 
orbital energy is within a large fragment of the tidally 
disrupted parent body. Surviving fragments may be re- 
sponsible for warps i nvoked in some di s k mod els, such as 
GD362 and GD56 (lJura et al.ll2009al l2007bD . Saturn's 
ring-moon Daphnis - 8 km in diameter and comparable 
to a small asteroid - causes 1 km vertical disturbances 
in the A ring; 100 times larger than the height of the 
unperturbed rings (Weiss et al. 2009). The more sub- 
stantial Mimas, lying just outside Saturn's main rings 
(and its Roche limit) with a mass nearly identical to the 
total mass of the rings, is responsible for the Cassini Di- 
vision between the A and B rings. A Mimas-analog of 
10^^ g could orbit a white dwarf beyond 2.0 Rq and not 
be tidally destroyed. Such a body with an inclined (rel- 
ative to the disk) or eccentric orbit, could easily excite 
warping or eccentricit y within a white dwarf dust ring 
(|Borderies et al.|[l984D . Dynamical modeling of tidal dis- 
ruption events may shed light on the frequency, size, and 
orbital distribution of residual fragments. 

5. CONCLUSIONS 

Three (possibly four) newly identified dust disks orbit- 
ing metal-polluted white dwarfs are discovered by Spitzer 
IRAC as mid-infrared excess emission over 3 — 8 fim. 
The fractional infrared (dust continuum) luminosity from 
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Table 8 

Metal-Polluted White Dwarf Target Parameters 



WD 


SpT 


M 




^cool 


[H/He] 


log < dhU/dt > 


log (Mz) 






(Mq) 


(K) 


(Gyr) 






(g) 


0047+190 


DAZ 


0.51 


16 600 


0.11 




9.0 




0106-328 


DAZ 


0.62 


15 700 


0.17 




9.3 




0307+077 


DAZ 


0.66 


10 200 


0.73 




8.9 




0842+231 


DBZ 


0.71 


18 600 


0.16 


-4.5=^ 


9.4 


21.7 


1011+570 


DBZ 


0.60 


18 000 


0.12 


-5.5=' 


8.7 


21.2 


1225-079 


DZA 


0.58 


10 500 


0.59 


-4.0 


9.2 


23.1 


1709+230 


DBAZ 


0.65 


18 500 


0.13 


-4.4 


8.4 


20.7 


2221-165 


DAZ 


0.73 


10100 


0.94 




8.9 




2229+235 


DAZ 


0.57 


18 600 


0.09 




9.3 




Note. — 


The seventh column lists 


time- averaged 


metal infall rates, 


calculated 



assuming a steady state between accretion and diffusion. The eighth column shows 
the mass of metals in the convective envelopes of the DBZ stars, based on the observed 
calcium abundances and assuming this re presents 1.6% of the total mass of heavy 
elements, as in the bulk Earth (Allegre et al.| 1993). 



^ Upper limit. 

HE 2221-1630 is relatively large and around 0.008. In 
contrast, HE 0106-3253 and HE 0307+0746 with their 
much smaller infrared luminosities are modeled by nar- 
row rings of orbiting dust. Even narrower disks at white 
dwarfs are possible, and perhaps likely given the neces- 
sity for ongoing heavy element pollution at the bulk of 
DAZ stars. An attenuated, flat disk as narrow as half 
an earth radius can harbor up to 10^^ g of material and 
be sufficient to reproduce the inferred 1) metal accre- 
tion rates of DAZ stars for at least 10'* yr, and 2) mass of 
heavy elements in the convective envelopes of DBZ stars. 
Yet such a tenuous ring of dust would not exhibit an in- 
frared excess at the 10% level for any inclination. The 
only evidence of accretion may be the photospheric pol- 
lution; other signatures such as X-ray emission seen in 
accreting binaries require much higher mass infall rates 
than inferred for metal-enriched white dwarfs. There- 
fore, elemental abundances for polluted stars without an 
infrared excess are still likely to represent the bulk com- 
position of accreted minor planets. 

PG 1225—079 is a tantalizing potential disk that emits 
only at 7.9 yum, but which may confirm a previously sus- 
pected class of dust ring with a large inner hole. While 
the excess emission from G166-58 is somewhat similar in 



temperature, a flat ring model predicts the inner hole at 
PG 1225—079 would be significantly larger and may rep- 
resent a disk that is near to being fully accreted; forth- 
coming data may test this scenario. 
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APPENDIX 

ACCRETION RATES FROM POYNTING-ROBERTSON DRAG ON DUST PARTICLES 



lJura et aD ()2007al ) demonstrate that dust captured near the Bondi-Hoyle radius (i.e. interstellar grains), and accreted 
under the influence of Poynting-Robertson (PR) drag, cannot account for the detected (and upper limit) infrared 
emission at DAZ white dwarfs. That is, assuming metal accretion rates sufficient to reproduce the photospheric 
metals, the predicted 24 fim emission from grains spiraling in ward under radiat i on drag is typical ly tens to hundreds 
of times higher than that measured at polluted white dwarfs (jFarihi et al.ll2009l : lJura et al.ll2007al) . 

Here, this argument is extended to consider PR drag on closely orbiting dust rings. To better understand the lifetime 
of such dust disks in the presence of accretion, consider the following simplistic but instructive model. A flat ring 
of dust extends outward from 0.2 i?Q, with thickness h = lOm (typical of Saturn's rings), and consisting of p = 2.5 
dust grains. T he ring orbits a compact star with L = 10~^ Lq, a value typical for a 16 000 K DA white dwarf 



cm 



Fontaine et al.|[2001|) . whose full starlight penetrates the in ner disk edge to a depth of one mean free path (r = 1). 
The PR lifetime of particles at the inner surface is given by ([Burns et al.lll979l ): 



Att c r'^pa 
3 LQp^ 



(Al) 



where r is the distance to the inner disk edge, a is the particle size, and Qpr is the radiative coupling coefficient. For 
particles around 1 ^m or larger, geometrical optics is valid and Qpr = 1. Applied to a disk of fixed volume and uniform 
density, the shell of mass exposed to the full radiation drag force is 
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Mpr = r;7r(r'2 - r^)h (A2) 

where 77 is the volume mass density of the disk, and r' — r is the depth of the exposed sheh mass (the mean free path 
in this case). For r' ^ |r' — r| One can make the approximation 



„'2 



r 



where / is the mean free path given by 



and the number density of particles, n is 



Hence the mean free path can be reduced to 



= (r - r){r' + r) !=i 2lr (A3) 



(A4) 



3^ 



(A5) 



(A6) 



and the exposed shell mass depends only on particle size, and is specifically independent of the disk density (or total 
disk mass in this case) 

Mp, = (A7) 

Because the PR lifetime of this shell is also proportional to particle size, the mass infall rate at the inner disk edge is 
essentially constant for the case of geometric optics 

M=^.^ (A8) 

and around 6000 gs^^ for the parameters given above. For silicate particles smaller than around 1 /im, geometric optics 
is not approp riate; radiative coupling becomes more effective until around 0.3 /im, then less effective at smaller radii 
(| Artvmowicz|[T988i) . The maximum implied mass infall rate due to PR drag calculated in this manner is 1.0 x 10''gs~^ 
and a few to several orders of magnitude shorter than the inferred metal accretion rates for DAZ white dwarfs. 

It should be noted that this model does not apply to disk particles excepting those at the innermost edge which are 
exposed to the full starlight. However, such particles should rapidly sublimate at some radius and are then no longer 
subject to radiation drag, further corroborating that PR forces cannot account for the necessary dust/gas infall rates. 
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